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^ ■ We show that translations and dilations of a p-adic wavelet coincides (up to the 

multiplication by some root of one) with a vector from the known basis of discrete 
p-adic wavelets. In this sense the continuous p-adic wavelet transform coincides with 
t*^ ■ the discrete p-adic wavelet transform. 

The p-adic multiresolution approximation is introduced and relation with the real 
multiresolution approximation is described. 
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O ' Relation of application of p-adic wavelets to spectral theory of p-adic pseudodiffer- 

. ential operators and the known results about sparsity of matrices of some real integral 

i operators in the bases of multiresolution wavelets is discussed. 
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^ : 1 Introduction 

Wavelet analysis (see e.g. [1], [2]) is an important method used in a wide variety of appli- 
cations from the theory of functions to signal analysis. For the standard wavelet analysis of 
functions of a real argument one can consider both the continuous wavelet transform and 
the discrete wavelet transform (expansion over wavelet bases). 

A basis of wavelets in spaces of complex valued functions of a p-adic argument was 
introduced and its relation to the spectral analysis of p-adic pseudodifferential operators 
was investigated in [3]. One can also consider, in analogy with the real case, a continuous 
p-adic wavelet transform, see |1] for a discussion. 

In the present paper we show that in the p-adic case, with the proper choice of a wavelet, 
an arbitrary translation and dilation gives (up to the multiplication by some root of one) 
some vector from the basis of discrete p-adic wavelets introduced in [3]. In this sense the 
continuous p-adic wavelet transform coincides with the discrete p-adic wavelet transform. 

This result looks strange from the point of view of real analysis. The observed behavior 
is a manifestation of the ultrametricity of the p-adic norm. In an ultrametric space it is not 
possible to leave a ball making steps smaller than the diameter of this ball. This implies the 
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existence of locally constant functions — functions which are constant on some vicinity for 
any point but not necessarily constant globally. 

When we apply the transformation of dilation or translation to a locally constant function 
of a p-adic argument, then, if the transformation is sufficiently small, the function will be 
invariant. This allows us to expect some properties of discreteness for continuous p-adic 
wavelet transform. In the present paper we show that making the right choice of wavelet it 
is possible to make continuous and discrete p-adic wavelet transforms to be equal. 

The next construction of the present paper is the multiresolution approximation of 
L^(Qp). Unlike in the real case, where translations by elements of the subgroup of integers 
in the real line were used, in the definition of the p-adic multiresolution approximation we 
use translations by elements of the factorgroup Qp/Zp. We show that the Bruhat-Schwartz 
space V{Qp) of p-adic test functions can be considered as a multiresolution approximation 
of L'^{Qp). In this sense the construction of multiresolution analysis (MRA) is a natural 
property of the p-adic analysis. 

These results, together with the simple relation of p-adic wavelets to spectral theory of 
p-adic pseudodifferential operators described in |3] (the real analogue of this relation we will 
discuss in Section 3), supports the point of view that the wavelet analysis is, in essence, a p- 
adic theory. In the real case we obtain a much more complicated wavelet analysis compared 
to the p-adic case. 

The structure of the present paper is as follows. 

In Section 2 we consider the continuous p-adic wavelet transform and show that it coin- 
cides with the discrete wavelet transform. 

In Section 3 we introduce the p-adic multiresolution approximation and discuss its rela- 
tion with the real multiresolution approximation. We also discuss relation of application of 
p-adic wavelets to spectral theory of p-adic pseudodifferential operators and the known re- 
sults about sparsity of matrices of some real integral operators in the bases of multiresolution 
wavelets. 

In Section 4, the Appendix, we put the material on p-adic analysis and p-adic wavelets, 
used in the present paper. 

2 Continuous p—adic wavelet transform 

For the notations used in the present section see the Appendix. 

Let us consider the main construction of the present paper — the continuous p-adic 
wavelet transform generated by the p-adic wavelet of the form 

ipix) = xip~^x)n{\x\p) 

which is equal to the product of a character x and the characteristic function Q of a. ball. 

In the next lemma we show that any vector from the basis {ip-ynj} of p-adic wavelets is 
a translation and dilation of the wavelet ip, and, vice versa, any translation and dilation of 
the wavelet ip belongs to the basis {ip'ynj} up to multiplication by some root of one. 

Equivalently, the basis of p-adic wavelets (multiplied by corresponding roots of one) 
can be considered as an orbit of action of the p-adic affine group in L'^{Qp) containing the 
function 
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Lemma 1 1) Any wavelet from the basis {ip-ynj} ofp-adic wavelets is a translation and 
dilation of the wavelet 

of the following form 

2) The translation and dilation of the p-adic wavelet ip is proportional to a vector from the 

basis {■07nj}' 

i/j'^'^x) = X {p~^ [(a|a|p)-^modp] {{\a\pb} - [|a|p6modp])) iJiog^\a\^,{\a\^b},{a\aM-^modp (1) 
Here a,b & Qp, a ^ 0. 
Proof 1) One has 



P ^x{P ip^jx ~ f^j))^{\p^jx — njlp) = p ^il^ip^jx — nj) = i/j- 



X] 



2) Let us compute 



X — b 



We have 



a — \a\ 



X — b 



fl {\\a\pX - \a\pb\p) = fl {\\a\pX - {\a\pb}\p) 



X \p ^ 



X — b 



X [P i-rd^lp^ ~ l^lp^) ^(Mp^ - {Hpb}\p) 



= (a|a|p)~^modp {\a\pX - \a\pb)'j Q {\\a\pX - {\a\pb}\p) = 

= X (p~^ [(a|a|p)"^modp ({|a|p6} - [\a\pb mod p])^ 

X {p~^ (a|a|p)~^modp] (\a\pX - {\a\pb})^ Q {\\a\pX - {\a\pb}\p) 
Taking into account the normahzation, we get 

= X {p~^ [(a|a|p)-^modp] {{\a\pb} - [\a\pb mod p])) ^iogp|a|p,{|a|p6},(a|a|p)-imodp 

i.e. the above translation and dilation of the wavelet ip coincides with a wavelet from the 
basis of p-adic wavelets times a root of one. This finishes the proof of the lemma. 

The next theorem shows that, in the p-adic case, the continuous wavelet theory in essence 
coincides with the discrete wavelet theory. 
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Theorem 2 1) The p-adic wavelet transform T, corresponding to the wavelet ip, takes 
the form of an expansion over the basis of p-adic wavelets ip^nj • 

(T/)(a,6) = |a|;^ / dxf{x)i^(^] = {rM) = 



= x{-P ^ [{a\a\p) ^modp] {{\a\pb} - [|a|p6modp])) f\og^\a\p,{\a\pb},{a\a\p)-^modp 

where f^nj is the coefficient of expansion of the function f over the basis of wavelets. 
2) For the continuous p-adic wavelet transform related to the p-adic wavelet ip the formula 
for the inverse wavelet transform takes the form of the expansion of the unit operator over 
the basis of p-adic wavelets 



7GZ,nGQp/Zp,j=l,...,p-l 



Proof The first statement foUows directly from Lemma [TJ 
Let us compute for the p-adic wavelet ip the constant 

Here we used ([T]) and the orthogonality of the wavelets in the basis, is the L^-norm. 
Then, using ([1]), we compute the integral 

dadb I, I, f dadb 

-— ' -Ji^ip ' \ - / -—^\'lpiog^\a\p,{\a\pb},{a\a\p)-^modp){W\ogp\a\^ - 
\^\p '' \"'\p 



H / ■r-rlV'logp|a|j„n,(a|a|p)-lmodp)(V'logj,|a|j„n,(a|a|p)-lmodp| = P ^ XI \'^inj){'^-^n]\ 

n<^Qp/Zp'' \^\P ■yeZ,neQp/Zp,j=l,...,p-l 



Therefore 

dadb 



This finishes the proof of the theorem. 
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7GZ,nGQp/Zpj=l,...,p-l 



3 Multiresolution approximation of L'^{Qp) and integral 
operators 

The following definition of a multiresolution approximation can be found in [2]. 

Definition 3 A multiresolution approximation of L'^{R) is a decreasing sequence Vy, 
^ & Z, of closed linear subspaces of L'^{R) with the following properties: 1) 

+ 00 +00 

Pi = {0}, y is dense in L'^iR); 

— oo — oo 
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2) for all f G L\R) and all e Z, 

/(■) EV,^ /(2-) e V,^^■, 

3) for all f e L^R) and all ne Z, 

f{.)eVo^fi--n)EVo; 

4) there exists a function g E Vq such that the sequence 

g{--n), neZ 

is a Riesz basis of the space Vq. 

It follows that there exists a function G Vq, such that the sequence 0(x — n), n E Z is 
an orthonormal basis of the space Vq. In the following we will speak only about ortho normal 
bases. 

We denote by the orthogonal complement to Vy in the space V^_i. Then L'^{R) is 
the orthogonal sum of the W^: 

L\R) = (B^ezW^ 

There exists a well known construction which allows to build wavelet bases from multireso- 
lution approximations, see [1] or [2] for the details. The space possesses the basis 
n G Z of multiresolution wavelets. 

Example 1 To define a multiresolution approximation it is sufficient to fix the function 
(p. Let us take this function to be equal to the characteristic function of the interval [0, 1]: 

(J^ix) = X[o,i]ix), xeR 

This choice of function (p is related to the Haar wavelets. It is easy to see that for the example 
under consideration all the properties of multiresolution approximation will be satisfied. 
Taking (j){x) = Xlo,i]{^) and n > in the definition above, we obtain the multiresolution 
approximation for the space L'^{R^) of quadratically integrable functions on the positive 
half-line. 

It is easy to find a p-adic analogue of the multiresolution approximation. The main 
difference is that in the p-adic case we should use translations by elements of the factorgroup 
Qp/Zp instead of translations by integers (we represent n G Qp/Zp by the numbers (jll), see 
the Appendix). 

Definition 4 A multiresolution approximation of L'^{Qp) is a decreasing sequence Vy, 
7 G of closed linear subspaces of L'^{Qp) with the following properties: 1) 

+ 0O +00 

Pi = {0}, U is dense in L'^{Qp); 

— oo —CO 

2) for all f G L\Qp) and all e Z, 

/(.) EV,^ fip-'-) G Vi; 
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3) for all f e L^iQp) and all n G Qp/Zp, 



f{.)eVo^f{--n)EVo; 
4) there exists a function G Vq such that the sequence 

0(--n), ne Qp/Zp 
is an orthonornial basis of the space Vq. 

We do not speak here about the Riesz bases since, as we will see, in the p-adic case we 
have a natural example of orthonornial basis satisfying the above properties. 

The space T>{Qp) of p-adic test functions possesses a natural filtration. Let us denote 
by V^{Qp) the space of locally constant compactly supported functions with the diameter 
of local constancy p"': 

\x - y\p <p^ ^ f{x) = f{y) 

We have the following theorem, which shows that the space L'^{Qp) possesses the natural 
multiresolution approximation by the space V{Qp) of test functions. Therefore the multires- 
olution approximation can be considered as an intrinsic property of spaces of functions of 
a p-adic argument. Moreover, the p-adic change of variable maps the multiresolution 
approximation of L'^{Qp) onto the multiresolution approximation of L'^{R+). 

Theorem 5 1) One has the filtration 

+ 00 

i^iQp) = U ^7(Qp); 

— oo 

2) The sequence V^{Qp), j E Z is a multiresolution approximation of L'^{Qp). For this 
multiresolution approximation the function 4>{x) can be taken to be equal to the characteristic 
function of the unit ball: 

(t>{x) = ^{\x\p) 

3) For p = 2 the map ([SP maps the multiresolution approximation {Vy} of LP'^Rj^) described 
in the Example 1 onto the multiresolution approximation {Vj{Qp)} of L'^{Qp). 

The corresponding isomorphism of spaces 'D^{Qp) and is given by the following for- 
mula: 

n{\p'^x - n\p) = X[o,p->]ipix) - p'^p{n)), x E Qp, nEQp/Zp. 

Proof The first statement of the theorem means that for any locally constant function 
with compact support there exists a nonzero infimum of diameters of local constancy. This 
follows from the compactness of the support. 

The proof of the second statements of the theorem is straightforward. The function (j) in 
the Definition H] can be taken equal to be to the characteristic function of the ball: 

= fi(|x|p) 

The translations Q{\x — n\p), n E Qp/Zp given by (jlj), provide an orthonormal basis in 
Vq^Qp). Dilations by the degrees of p give the spaces T>^{Qp). 
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The third statement follows from Lemma [7] and the above choice of the function 0. The 
space V^{Qp) possesses the basis — n|p)}, x G Qp, n G Qp/Zp. The set of functions 

{X[Q,p-i]{.^ ~ '^)} 1 ""^ ^ ^0, X G R+ constitutes a basis in V^. Since the map p is a one-to-one 
correspondence between Qp/Zp and A^o^ this finishes the proof of the theorem. 

Note that, contrary to the real case, where characteristic functions and Haar wavelets 
are discontinuous, the multiresolution approximation of L^{Qp) considered in the above 
theorem consists of spaces of continuous functions. After the p-adic change of variable the 
construction which was not regular (was discontinuous) in the real case, becomes highly 
regular in the p-adic case. 

Discuss now application of wavelets to investigation of integral operators. It was found 
[5] that for some families of integral operators, say the Calderon-Zygmund and the pseu- 
dodifferential operators, matrices of these operators will be close to diagonal in wavelet 
bases for which wavelets has many momenta equal to zero. Consider the multiresolution 
approximation for L^{R) for which the corresponding wavelet has M vanishing moments: 

J tp{x)x"'dx = 0, m = 0,l,...M-l. 

Then [5j we have the following property of decay for matrix elements of some integral 
operators 

+ m + < 1 + |f "^|M^^ 

for all 

\i-l\> 2M 

Here the matrix elements a]i, f3Ji, are taken in the wavelet bases in the spaces Vj and 
W-y of the multiresolution approximation, the index 7 is related to the scale of the wavelets 
and the indices i, I related to space localization (i.e. these are translation indices) of the 
wavelets. 

This means that matrices of the corresponding integral operators in the wavelet bases are 
sparse and close to diagonal matrices. This result has a variety of applications, for example 
to numerical algorithms of fast multiplication of integral operators in wavelet bases. 

Discuss now the p-adic analogue of this result. It is known [3], see the Theorem [6] in 
the Appendix, that p-adic pseudodifferential operators are diagonal in the basis of p-adic 
wavelets. Moreover, in the p-adic case we do not need any conditions of vanishing of higher 
moments for the wavelets. Therefore in the p-adic case we have much stronger results about 
the relation of wavelet bases and integral operators. The diagonality of matrices of integral 
operators which was approximate for the real case becomes exact for the p-adic case, and, 
moreover, we are able to compute the spectra of the corresponding operators. 

4 Appendix 

Let us recall some constructions of p-adic analysis, see [6] for details and, e.g. [7], [8] for 
applications. The field Qp of p-adic numbers is the completion of Q with respect to the 
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p-adic norm | ■ \p, defined as follows. For any rational number we consider the representation 

X = p' — 
n 

where p is a prime, 7 is an integer, p, m, n are mutually prime numbers, n 7^ 0. The p-adic 
norm is defined as follows 

\x\p =p~^ 

A j9-adic number can be uniquely represented by the series 

00 

a; = ^ Xjpr' , = 0, . . . , p — 1 

i=7 

which converges in the p-adic norm. 

A complex valued character x : Qp ^ C of a p-adic argument (where x has the form of 
the above series) is defined by 

X{x) = exp ^Tci J2 ^jP'^ , X{x + y) = x{x)xiy) 

When 7 above is nonnegative the sum above is equal to zero and the character is equal to 
one. 

The character x is a locally constant function. The function / is locally constant if for 
any x there exists e for which \/y: \x — y\p < e we have f{x) = f{y). 

Another example of a locally constant function is characteristic function of a p-adic ball. 
We denote fl{x) the characteristic function of the interval [0, 1]. The characteristic function 
of the p-adic ball of radius 1 with center in has the form: 



n{\x\ 



\ 1. 


\x[ 


10. 


\x[ 



p < 1 
p > 1 



Let us note that, unlike in the real case, the characteristic function of a p-adic ball is 
continuous (the same holds for an arbitrary locally constant function). 

The Bruhat-Schwartz space T>{Qp) of p-adic test functions is the linear space of locally 
constant complex valued functions with compact support. Any function in T>{Qp) is a (finite) 
linear combination of characteristic functions of balls. 

The Vladimirov operator of p-adic fractional differentiation has the following form 

D f{x) = — / -—dfi{y), a>0 

Here /z is the Haar measure on p-adic field (for which the measure of a ball is equal to the 
diameter of this ball). The p-adic gamma function has the form 

The domain of contains the space T>{Qp). 

The Vladimirov operator generates a natural random walk on the p-adic field considered 
in [7j . For more general considerations of p-adic dynamics see [8] . 

The basis of discrete p-adic wavelets and its relations to the spectral theory of p-adic 
pseudodifferential operators is described in the following theorem [3]: 
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Theorem 6 1) The set of functions {ip^^/nj}'- 

ip^njix) = p'2^(p^-^j(x -p^'^n))n{\p'^x - n\p), 

xeQp, 7 e Z, neQp/Zp, j = l,...,p-l (2) 

is an orthonormal basis in L'^{Qp). 

2) This basis consists of eigenvectors of the Vladiniirov operator D": 

/^"^-yn.- = p"^'-^V7ni (3) 

here n G Qp/ Zp in (0) are represented by the numbers 

n = ^nip\ ni = 0,...,p-l (4) 

l=f3 

It is possible to compute, using the basis of p-adic wavelets, spectra of a very wide 
class of p-adic pseudodifferential operators [9] . Multidimensional p-adic wavelets in relation 
to spectral theory of p-adic pseudodifferential operators were discussed in [TU] . Wavelets 
on some family of abelian locally compact groups were considered in [TT]. There exists a 
generalization of theory of p-adic wavelets and pseudodifferential operators onto general 
locally compact ultrametric spaces [12], [13], [Hj. 

There exists a natural map of p-adic numbers onto the positive half-line i?+ = [0, oo), 
which, for p = 2, maps the basis of 2-adic wavelets onto the basis of Haar wavelets on the 
half-line. This exhibits a natural relation between real and j9-adic wavelet theories. 

Let us discuss this relation. The wavelet basis in L'^{R) contains the functions 

= 2"i^(2"^x-n), 7GZ, uEZ 
where \l/ is some integrable mean zero function. The simplest example is the Haar wavelet 

equal to the difference of two characteristic functions. 

Define the p-adic change of variable (or the Monna map) as follows: 

p:Qp^R+ 

OO OO 

p : ^Xip* ^Xip"'~\ Xi = 0,...,p- 1, 'J e Z (5) 

i=7 i=7 

This map is one-to-one almost everywhere and conserves the measure: 

/i(S) = L{p{S)) 

where S" is a measurable set, /i is the p-adic Haar measure, L is the Lebesgue measure on 
the half-line. The map p satisfies the Holder inequality: 

\p{x) - p{y)\ < \x-y\p 

The p-adic change of variable is a one-to-one map between the group Qp/Zp and the set 
Nq of natural numbers including zero. Here the elements of Qp/Zp are given by the rational 
numbers (jlj). 

The following lemma and theorem were proven in [5]. 
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Lemma 7 For n G Qp/Zp and m, k & Z , k > m the map p satisfies the conditions 

p:f^n + p'^Zp p-™p(n) + [0, p-'^] (6) 

P : QpMv'^n + p'Zp} i?+\{j9->H + (0,^-^=)}, n ^ (7) 
p:QA{P%}^^A[0,p-') 



Theorem 8 For p = 2 the map 1^ provides a one to one correspondence between the 
Haar basis in L'^{R^) and the basis of p-adic wavelets in LP'{Qp): 

'^-tp(n){p{x)) =i)^nl{x) (8) 

We understand formula ([S]) in the sense (since the map p is not a one-to-one map on 
the set of measure zero). 

One can consider the continuous wavelet transform on the field of p-adic numbers [1]. 
The affine group on the field of p-adic numbers in the unitary representation G in L'^[Qp) 
acts by translations and dilations 

G{aMf{x) = r\x) = \a\7'f[^^ 

where a,b E Qp and a 7^ 0. 

Assume we have a wavelet on Qp — i.e. a mean zero complex valued function in 
L^iQp) r\L'^{Qp)- The wavelet transform Tf of the function / G L'^{Qp) is defined by 

(T/)(a,6) = |a|;^ / dxf{x)^{^\ = {^^'\ f) 



where (■, ■) is the scalar product in L'^{Qp). 
The inverse wavelet transform has the form 



fix) = C^' I ^{^''^ m'^^'ix) = C^' I ^iTf){a, b)^'^'\x) (9) 

I ^ I J) I ^ 



Here 

dadb 



l^^P J \a\ 
■ I II I 1^ 

and is the L^-norm. 

The formula for the inverse wavelet transform (Q has the equivalent form of an expansion 
of the unit 



Here we use the standard notations from quantum mechanics: IvEf"'*) denotes the element 
of the Hilbert space L^(Qp) given by the function vlf^'*, is the canonically conjugate 

linear bounded functional on L?{Qp), and |^°''')(vl/'^'*| is the orthogonal rank one projection 
onto ^"'^ G L'^{Qp). 
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